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Abstract. We report on two photometric monitoring campaigns of Very Low Mass (VLM) objects in the young open cluster 
around crOrionis. Our targets were pre-selected with multi-filter photometry in a field of 0.36sqdeg. For 23 of these objects, 
spanning a mass range from 0.03 to 0.7 Mq, we detect periodic variability. Of these, 16 exhibit low-level variability, with 
amplitudes of less than 0.2 mag in the I-band, which is mostly well-approximated by a sine wave. These periodicities are 
probably caused by photospheric spots co-rotating with the objects. In contrast, the remaining variable targets show high- 
level variability with amplitudes ranging from 0.25 to 1.1 mag, consisting of a periodic light variation onto which short-term 
fluctuations are superimposed. This variability pattern is very similar to the photometric behaviour of solar-mass, classical T 
Tauri stars. Low-resolution spectra of a few of these objects reveal strong Ha and Ca-triplet emission, indicative of ongoing 
accretion processes. This suggests that 5-7% of our targets still possess a circumstellar disk. In combination with previous 
results for younger objects, this translates into a disk lifetime of 3-4 Myr, significantly shorter than for solar mass stars. The 
highly variable objects rotate on average slower than the low-amplitude variables, which is expected in terms of a disk-locking 
scenario. There is a trend towards faster rotation with decreasing mass, which might be caused by shortening of the disk 
lifetimes or attenuation of magnetic fields. 
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1. Introduction 

Open clusters are ideal environments to study stellar proper- 
ties and evolution, because they contain a homogeneous pop- 
ulation of objects with known distance, metallicity, and age. 
Recent deep surveys have unveiled the population of several 
open clusters far down into the substellar regime. Examples 



are the surveys by Moraux et al. (2003 1, Pinfield et al. (2000 1 
and Zapatero Osorio et al. (1999 1997 1 in the Pleiades, and 
the work of Barrado y Navascues et al. in IC2391 (2001 1 



and a Per (2002 1. In two cases, namely in the clusters around 
a Ori (Zapatero Osorio et al. 120001 and the Orion Trapezium 
Cluster (Lucas & Roche 2000 1, even isolated planetary mass 
objects were found. Thanks to these surveys, large samples of 
Very Low Mass (VLM) objects are known today, among which 
we conveniently will subsume all objects with masses below 
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0.4 Mq, including very low mass stars, brown dwarfs and free- 
floating planetary mass objects. All objects with masses below 
this 0.4 Mq limit are thought to be fully convective (Chabrier 
& Baraffe 2000 1, giving a physical motivation for this defini- 
tion of VLM objects. 

Wide-field photometric monitoring is a powerful tool to in- 
vestigate properties of large object samples: If an object ex- 
hibits asymmetrically distributed surface features, e.g. magnet- 
ically induced star spots, its flux will be modulated with the 
rotation period. Hence, period search in the lightcurve allows 
the determination of the (projection-free) rotation period. The 
amplitude of the periodicity, in turn, depends on the properties 
of the star spots, allowing conclusions about surface activity 
processes. Pointing to dense open cluster fields, such monitor- 
ing campaigns become very efficient, as one can register many 
objects contemporaneously within one field of view. 

Compared to solar-mass stars, there is a significant lack of 
known rotation periods for VLM objects. Photometric moni- 
toring studies delivered a small number of periods for evolved 
VLM stars and Brown Dwarfs (Bailer- Jones & Mundt IT9991 
I200T1 Ti nney & Tolle y [19991 Martin et al. I200T1 Clarke et 
al. 120021 Gelino et al. 2002 1. Only three periods are known 
for VLM objects with ages between 50 and 100 Myr (Martin 
& Zapatero Osorio [T997I Terndrup et al. 1999). Two periods 
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for VLM members of a Ori were published by Bailer- Jones 
& Mundt (1200 1> . Recently, Joergens et al. (2003) report on 5 
measured periods for Brown Dwarfs and VLM stars in the very 
young Chal star forming region. Our own study in the young 
open cluster IC4665 (age 36 Myr) delivered several rotation pe- 
riods for VLM objects (Eisloffel & Scholz 120021 

Observations of very young VLM objects are of special 
interest, because they can deliver clues about their formation 
process. Recent results suggest that objects in the substellar 
regime form similar to stars. Several authors detected typical T 
Tauri star phenomena on VLM objects, e.g. outflow processes 
(Fernandez & Comeron 2001 1 and mid-infrared excesses at- 
tributed to the presence of a circumstellar disk (Natta & Testi 
2001 Apai et al. l2002l Jayawardhana et al. 120031 . Moreover, 
Lopez Martf et al. (2004» and Liu et al. {2003 1 demonstrate a 
correlation between Ha emission and mid-infrared excess for 
Brown Dwarfs down to 0.02 M Q , which they interpret as an 
indication for ongoing accretion. 

Solar-mass T Tauri stars show various types of photomet- 
ric variability. Active accretion processes often manifest them- 
selves by large amplitude variations showing several possible 
periods (e.g., Fernandez & Eiroa 1996 Bouvier et al. 1995 
Herbst et al. 120001 . If VLM objects undergo a T Tauri phase 
as well, they should exhibit a similar photometric behaviour. 
Hence, photometric monitoring can deliver an independent 
contribution to the ongoing debate about VLM object forma- 
tion. 

The a Ori cluster is an appropriate target for such a study, 
because it is relatively nearby (350pc, Bejar et al. 1999 1, has 
negligible extinction (£b-v = 0.05, Bejar et al. I1999i and an 
age of 3 Myr (Zapatero Osorio et al. 2002 1. The extended work 
by Bejar et al. and Zapatero Osorio et al. revealed a rich VLM 
population. We monitored this cluster in two photometric time 
series. Complementary observations, presented in Sect.|3 iden- 
tified the cluster members in the time series field. We report on 
the monitoring campaigns in Sect.|3]and the time series analy- 
sis in Sect.|4] The following sections describe the results of out 
lightcurve analysis: We first establish two origins for the ob- 
served variability (Sect.|5Jl. In Sect.|5J we then concentrate on 
the investigation of the rotation periods. SectionQcontains the 
discussion of low-resolution spectra for highly variable objects. 
Finally, in Sect.|8] we present our conclusions. 

2. Selection of VLM objects in the a Ori cluster 

We searched for VLM members of the erOri cluster using 
the 2-m Schmidt camera at the Thiiringer Landessternwarte 
Tautenburg (TLS) and the 1.23-m telescope on Calar Alto 
(CA) in combination with infrared photometry from 2MASS 1 . 
Figured shows the 36' x 36' TLS field and the 17' x 17' CA 
field in the north of a Ori. These survey images were obtained 
as part of the time series campaigns (see next section). The CA 
field was observed in the R- and I-band with exposure times of 
30 and 600 sec, respectively. For the TLS field, we used a 600- 
sec-exposure in the I-band. These images were reduced with 

1 Catalogue available under http://www.ipac.caltech.edu/2mass 
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Fig. 1. Time series fields in the a Ori cluster: The TLS field is 
indicated with solid, the CA field with dashed lines. The posi- 
tions of the cluster member candidates (see Table 7) are over- 
plotted, a Ori is the bright star near the bottom of the field. 

the standard routines within the ccdred package of IRAF 2 , in- 
cluding overscan and bias subtraction, as well as flatfield cor- 
rection. The I-band fringe pattern caused by night sky emission 
was corrected with a fringe mask constructed from blank field 
images. The sigma-clipped median of these frames ('superflat') 
was smoothed and subtracted from the original superflat, deliv- 
ering the required fringe mask. Because the fringe mask con- 
tains only an average fringe pattern, this mask must be scaled to 
properly subtract the fringes in the individual images. The scal- 
ing factor was determined as the ratio of the fringe amplitude 
in the fringe mask and in the respective image. 

The pixel positions of the objects in the I-band images were 
determined with the SExtractor software (Bertin & Arnouts 
1996 1. The respective positions in the R-band image were 
found by a linear transformation determined with several bright 
stars whose pixel positions were manually determined in R 
and I. The pixel positions were transformed to sky coordinates 
with the known sky coordinates of unsaturated HST guide stars 
(Morrison et al. 1200 1 1 in our field using the ccmap routine of 
IRAF. The coordinate precision is zbl'.'O. 

Instrumental magnitudes for all objects in the TLS and 
CA catalogues were determined by PSF fitting (TLS) or aper- 
ture photometry (CA) using the daophot package within IRAF 
(Stetson 1987). The PSF fitting photometry was considered to 
be more appropriate for the TLS data, because the poor seeing 
turned the cluster into a 'crowded field' . Near-infrared photom- 
etry in J, H and K for both fields was obtained from the 2MASS 
database. 

2 IRAF is distributed by National Optical Astronomy 
Observatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the National 
Science Foundation. 
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Fig. 2. (1,1- J) colour magnitude diagram for the CA (upper 
panel) and the TLS field (lower panel). All probable cluster 
members are marked as larger dots. Error bars indicate typi- 
cal photometry errors for the candidates. The separation line 
between cluster member candidates and field objects is shown 
as a straight line. The position of the 3 Myr isochrone from 
Baraffe et al. (1998 1 is indicated by the right-hand line. The 
arrow shows the reddening vector for Ay — 1 mag. 

In order to calibrate the CA images, we observed two stan- 
dard fields (Landolt 1 1992> at an airmass comparable with that 
of the a Ori field. I-band zero-point and colour coefficient were 
derived fitting the instrumental magnitudes and the catalogue 
magnitudes with the following relation (I - Landolt magnitude; 
i,r - instrumental magnitudes): 

/ = ica + ZP + C{r CA - ica) (1) 

With the resulting parameters (ZP = 3.45, C = 0.171), the 
a Ori magnitudes were converted to the Landolt system. The 
high colour coefficient points at significant differences between 
the CA and the Landolt I-band. Since the Landolt standards are 
mainly stars with R-I<1.0, one has to be cautious to apply this 
colour correction to very red targets. For such objects, we may 
overestimate the I-band flux. 

The CA field completely overlaps with the TLS field. 
Therefore, the TLS magnitudes were calibrated using objects 
contained in both catalogues. We noted a significant colour de- 
pendency of the zero-point offset between CA and TLS I-band, 



therefore we applied a colour correction using the 2MASS J- 
band magnitudes. The following transformation was used to 
convert TLS magnitudes into the Landolt system. 

I = iTLS + ZP + C{i TLS - J) (2) 

After this procedure, I-band magnitudes from both telescopes 
are available in the same photometric system. 

Cluster member candidates were selected from (I,TJ) 
colour magnitude diagrams (Fig.|2}. In both diagrams, a cu- 
mulation of objects to the right of the field stars is clearly visi- 
ble. To determine the position of this cumulation precisely, we 
divided the colour magnitude diagrams (CMD) in horizontal 
bins of Al = 1 mag. The histograms of the I-J colour of these 
bins typically show a broad maximum on the left and a second 
smaller maximum on the right side. This second peak indicates 
the position of the cluster isochrone. We registered the I-J value 
of this peak for each bin in both diagrams and fitted these val- 
ues linearly. This fit, the empirical isochrone, shifted 0.3 mag to 
the left to account for photometric errors and the uncertainty of 
the derived isochrone, was used as dividing line between clus- 
ter member candidates and field stars (see Fig.|2] straight line). 
All objects to the right side of this separation line are potential 
low-mass cluster members. Some of these objects clearly ap- 
pear reddened in the CMD. Since extinction towards a Ori is 
negligible (Bejar et al. 1999), this is a sign for intrinsic redden- 
ing, as found for objects with similar masses in the same cluster 
(Oliveira et al. 2002 1. We selected 55 candidates from the CA 
and 112 from the TLS diagram. The TLS candidates contain 28 
objects also identified with the CA photometry. The remaining 
CA candidates are either too bright for the TLS photometry or 
not resolved in the TLS image. 

In Fig.|2l we also indicate the position of the Baraffe et al. 
i 19981 1 isochrone for the a Ori cluster. The I-J colours of our 
candidates are significantly smaller than those from the evolu- 
tionary track for the appropriate age. This difference is possi- 
bly caused by insufficient colour correction in the I-band (see 
above) as well as shortcomings in the modeled BD atmospheres 
and colours. However, it is valid to compare the 2MASS pho- 
tometry with the Baraffe et al. isochrone. Systematic mis- 
matches between the 2MASS photometric system and the CIT 
system used for the model isochrones are safely below 3% 
(Carpenter 2001 ) and thus negligible for our purposes. 

We used the R-band photometry from CA to verify the sta- 
tus of the CA candidates. With one exception (which was re- 
jected), all CA candidates with detection in the R-band image 
are situated on the red side of the (LR-I) CMD. As a further 
test, we compared the (J-H) colours of our candidates with the 
models of Baraffe et al. ( 1998 1. The (J-H) values were preferred 
over the (J-K) colour because the latter might be influenced by 
excess radiation from a circumstellar disk, as already found for 
VLM objects of similar age (Muench et al. 120011 Oliveira et 
al. I20021 I. Figure0 shows the (J,J-H) CMD for all our candi- 
dates and the Baraffe et al. isochrones for the a Ori cluster. Our 
VLM candidates fall around the 3 Myr isochrone (solid line). 
The 1 Myr (dotted line) and 5 Myr (dashed line) isochrones are 
indistinguishable from the 3 Myr track for J > 12 mag, but 
show clear deviations for brighter objects. We excluded one 
CA and two TLS candidates which are clear outliers and thus 
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Fig. 3. (J,J-H) colour magnitude diagram for all candidates 
(plotted as filled squares). Error bars indicate typical photom- 
etry errors for the candidates. The 3Myr isochrone for cOri 
(Baraffe et al. 119981 is shown as a solid line. For bright ob- 
jects, the isochrones for 1 Myr (dotted line) and 5 Myr (dashed 
line) differ significantly from the 3 Myr track and are therefore 
also shown. Objects whose J-H colour deviates clearly from 
the cluster isochrone are marked with crosses and excluded. 
The positions of 1000 arbitrarily selected field stars are shown 
as small dots. 



probably non-members of the cluster. It should be noted that 
the (J,J-H) CMD is not particularly well suited to discriminate 
between cluster members and field stars, as seen in Fig.[3]from 
the 2MASS data for arbitrarily selected objects in the region 
around erOri. In contrast to the (1,1- J) CMD, there is no clear 
separation of field stars and cluster members. Therefore, the 
(J,J-H) CMD will only reveal obvious false detections and does 
not signficantly reduce the contamination rate of our sample 
(see below). 

The remaining candidate list comprises 52 CA and 110 
TLS objects, including 27 double detections. Coordinates and 
photometry for these 135 objects are listed in Table 7. Masses 
for these candidates were estimated by fitting the 3 Myr 
isochrone of Baraffe et al. for the J- and H-band with a 
low degree polynomial and applying the fit to the J- and H- 
magnitudes. This gives us two mass estimates, one based on J-, 
the other on H-band photometry. We found both values to be 
very similar, and use their average as final mass (see Table 7). 
The theoretical isochrone covers the mass range from 0.02 to 
1 .4 Mq . The mass estimate for fainter or brighter objects is thus 
not possible with these models. For these (very few) targets, we 
only give an upper or lower limit for the mass. 

Our candidate list may be contaminated by fore- and back- 
ground objects. However, the spectroscopic follow-up observa- 
tions of Bejar et al. (2001 1 for an object sample selected with 
similar photometric criteria suggest that contamination rates 
are low. A 7 arcmin high stripe in the southern part of our field 
is also covered by the survey of Bejar et al., although they are 
only sensitive to objects with / > 15. In the overlapping re- 
gion, we found 18 cluster member candidates in this magni- 
tude range. Of these candidates, 1 1 (61%) are also identified by 



Table 1. Cluster member candidates which were previously 
identified by Bejar et al. (2001 ) as possible a Ori members. We 
also list the I-band magnitudes given by Bejar et al. (2001} as 
well as the difference to our I-band photometry. 



No. 


Name 


Io (mag) 


Al 


4 


SUri JID3813. 1-Uzz41l) (SUnlj) 


10.41 


U.68 


9 


SOri J053817. 1-022228 (SOri9) 


15.81 


0.56 


22 


SOri J053829.5-022517 (SOri 29) 


17.23 


0.66 


33 


SOri J053849.2-022358 


16.81 


0.98 


49 


SOri J053915. 1-022152 (SOri 38) 


17.64 


0.10 


73 


SOri J053848.0-022854 (SOri 15) 


16.79 


0.83 


81 


SOri J053907.4-022908 (SOri 20) 


17.32 


0.90 


83 


SOri J053907.9-022848 (SOri 8) 


15.74 


0.37 


86 


SOri J053911.7-022741 (SOri 1) 


15.08 


0.36 


106 


SOri J053944.4-022445 (SOri 10) 


16.08 


0.34 


108 


SOri J053948. 1-022914 


18.92 


0.82 



Bejar et al. (2001 1. Assuming the completeness of their survey, 
this gives a first rough estimate of up to 40% for the contami- 
nation rate. Table^contains our identification numbers as well 
as object names and I-band photometry of Bejar et al. (|2001 1 
for common candidates. Our I-band magnitudes for these ob- 
jects are systematically smaller than those given in Bejar et al. 
( 120011 ). The offset between both photometric systems (given in 
the last column of TableQ increases with the I-J colour, which 
may result from insufficient colour correction in the absolute 
photometry. Since the objects in Tabled cover only a small 
range in I-J (8 of 11 objects have 0.4 < 7 - J < 0.7), it is 
not possible to derive a colour correction that would bring all 
objects onto the same system. 

An alternative estimate for the number of contaminating 
field stars can be obtained with the Besancon Galaxy model 3 
(Robin & Creze 1986). The simulation delivers star counts de- 
pending on Galactic coordinates, colour and brightness. A sim- 
ulated (I,R-I) CMD for a 1 sq field centred on a Ori contains 90 
objects in a 0.4 mag wide stripe around the cluster isochrone of 
Baraffe et al. ( 1998). Scaling to the appropriate field sizes, we 
calculate that 7 of the 52 CA candidates (13%) and 32 of the 
110 TLS (29%) candidates are field stars, i.e. a significantly 
lower rate than in our first estimate. Since interstellar extinc- 
tion variations might hamper the CMD simulation, we assume 
conservatively that about 30% of our candidates may not be 
members of the erOri cluster. However, the definite decision 
about the cluster membership of our candidates must be post- 
poned until we have obtained follow-up spectroscopy for all of 
them. Our preliminary cluster member list in Table 7 will be 
used in the following as target list for the variability study. 

3. Time Series Observations 

3. 1. Observing campaigns and data reduction 

The region with our cluster member candidates was monitored 
in two I-band time series campaigns, with the TLS Schmidt 

3 These simulations are available online under 

http://www.obs-besancon.fr/www/modele/modele_ang.html 
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telescope and with the 1.23-m telescope on Calar Alto (CA). 
Details of these runs are given in Table|2] During the CA cam- 
paign, we took alternating short (30 sec) and long (600 sec) ex- 
posures to extend the dynamical range of the time series pho- 
tometry. In both campaigns, we pointed to the same position in 
the sky to within w 10" in all monitoring nights. 

Table 2. Monitoring campaigns 





TLS Schmidt 


CA1.23 


aperture 


1.34 m 


1.23 m 


CCD size 


2kx2k 


2kx2k 


pixel scale 


l'.'26/pix 


0"5/pix 


observing nights 


16-26 Jan 2001 


16-27 Dec 2001 


mean seeing 


2'.'5 


1'.'8 


CCD images 


77 


78/78* 



* The CA time series consists of 78 short and 78 long exposures, 
which were taken alternately. 

Both observing runs span at least ten days. The distribution 
of the datapoints over the observing run, however, is irregular 
because of variable weather conditions. The number of data- 
points per night varies from to 27. Figure|4] shows the data- 
point distribution for both runs. 

Image reduction and photometry of both campaigns was 
done following the recipes explained in the previous section. 
An object catalogue was created for the deepest time series im- 
age of both campaigns. For each time series image, we deter- 
mined the spatial offset between the object catalogue and the 
positions in the actual image. Applying this offset to the object 
catalogue and centering these new positions on the respective 
image generates an object catalogue for each image. We used a 
fixed sample of PSF stars for each night of the TLS time series. 
After the photometry procedure, we obtained lightcurves for all 
objects in the respective object catalogue. 

3.2. Relative calibration 

The lightcurves of both campaigns were calibrated using an 
average lightcurve from a set of non-variable reference stars 
in the observed field. This is necessary to correct for changing 
atmospheric conditions as well as the variable airmass over the 
time series. In the following, we describe the selection process 
for the reference stars. 

First, a sample of stars was selected, whose photometry 
errors were below 0.1 mag in all images. To reject variable 
objects from this initial reference star list, we examined the 
variability of each star with respect to the average lightcurve 
of all other stars. We used the following procedure based on 



the routine by Allain ( 1994 1. The number of a reference star 
is indicated with i — 1 . . . Nr, the number of an image with 
j = l...N B . 



1 . averaging the instrumental magnitudes for each star 

, n b 
m~= Tr-Y^fe) 

3=1 



(3) 




Days from 16 Jan, UT12:00 




Days from 16 Dec, UT 12:00 

Fig. 4. Datapoint distribution for the TLS (upper panel) and the 
CA run (lower panel). Plotted is the non-integer fraction of the 
observing times against the observing times. 

2. subtraction of this average from every value of the time se- 
ries for every star 



(4) 



3. calculation of average m Q - and standard deviation <Tj for the 



differences from Eq. (4) 



, N R 



\ 



Nr 



Nr 
i=l 



(mHtj 



m 0)2 



(5) 



(6) 



4. rejection of obvious 'bad' images by examination of dj 
(N' B is the number of the remaining images) 

5. repeated execution of steps 1-3 without the rejected images 

6. calculation of a quality number testi, which allows for 
every star an examination of its intrinsic variability com- 
pared to all other reference stars; if \rn®(tj) 



testi j = 1> elsewhere 



N' B 

testi = test 



(7) 
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After the described procedure, the stars with test j < 3 were 
selected as reference stars. We obtained « 170 reference stars 
for the TLS and « 80 for the CA campaign. For all images, the 
average brightness of the reference stars was determined. 

N re{ ^ 

This mean lightcurve was subtracted from all time series. Thus, 
we obtained extinction corrected relative magnitudes m rel (tj) 
from the instrumental magnitudes for all objects. m(tj). 

m rel {tj) = m(tj) - m r ^{t 3 ) (9) 



with many bright stars lead to many saturated pixels which 
are useless for difference image analysis and decrease the 
precision of the fitting process. 
- Every bright star causes a so-called Schmidt ghost in the 
image, disturbing the fitting process. However, if the ghosts 
were masked so that the fit can be done without them, we 
would loose a considerable fraction of the targets. 

For all these reasons, we returned to 'classical' PSF fitting pho- 
tometry for the TLS images. High background caused by light 
reflection from snow, bad seeing and variable atmospheric con- 
ditions limited the precision of the TLS lightcurves. Here, we 
reach 0.015 mag for the brightest stars. 



3.3. Optimal Image Subtraction 

For CCD images from the 1.23-m telescope on the Calar Alto, 
a dedicated image analysis package has been developed by 
the Wendelstein Calar Alto Pixellensing Project (WeCAPP, 
Riffeser et al. 2001) team. This software is based on 'Optimal 
Image Subtraction' (OIS, Alard & Lupton [T9"98l . After astro- 
metric and photometric alignment, it performs an image con- 
volution with respect to a reference frame (Gossl & Riffeser 
2002 1. The reference image is obtained by coaddition of the 
time series images with the best seeing. For the long exposure 
time series from CA, we used this pipeline to check the quality 
of our lightcurves from our already described 'classical' pho- 
tometry. Because of bad image quality, six images were not 
used for this procedure. 

The basic idea of difference image analysis algorithms like 
OIS is to measure only the variability and not the constant por- 
tion of the star's flux. The result of the pipeline are frames 
which contain only variable sources. On these frames, we per- 
formed aperture photometry for the previously determined ob- 
ject catalogue. The resulting fluxes were transformed to rela- 
tive fluxes by division by the flux of the respective object in the 
reference image. Relative fluxes were transformed to relative 
magnitudes with rn re i = — 2.51og 10 (l + / rc i). 

To compare the results from OIS with those from daophot, 
we determined the photometric precision of both measure- 
ments: After excluding 3a outliers, we calculated average and 
rms for all lightcurves. The lightcurves of the brightest objects 
scatter with 4 mmag in the OIS data, 3 mmag less than with 
classical photometry. Throughout the usable magnitude range, 
OIS improves the precision by several mmag. Therefore, we 
used the lightcurves from the WeCAPP software for the fur- 
ther analysis. The short exposure CA images, however, cannot 
be reduced properly with the OIS pipeline, since they do not 
contain enough stars for the fitting process. 

Three obstacles prevented a successful application of OIS 
to the TLS data: 

- The pixel scale of the TLS Schmidt camera is roughly two 
times that of the CCD at the 1.23-m telescope on CA. This 
leads to an undersampled PSF, which makes exact image 
centering and image folding difficult. 

- The sky background in the TLS images was in most cases 
> 25000 ADU, because of snow cover around the obser- 
vatory. Therefore, moderately crowded open cluster fields 



4. Time Series Analysis 

4. 1. Generic variability test 

The candidate lightcurves of both campaigns were analysed 
with a unique procedure. In a first step, we examined all 
lightcurves visually to register obvious signs of rapid variabil- 
ity. In particular, we looked for sudden brightness eruptions 
(e.g. flares), eclipses or totally irregular variability. This search 
was not successful. It turned out that some TLS objects are 
either too bright or too faint to be detected in most of the 
images; time series analysis is not possible for these targets. 
Subsequently, the lightcurves were filtered: In three iterations, 
we excluded 3er outliers, but only if they were framed by other 
datapoints to avoid the exclusion of intrinsic short-term vari- 
ability. The root mean squares of the filtered lightcurves are 
listed in Table 7 and plotted in Figs.[5]and[6] The solid lines in 
these Figures indicate our photometric precision determined by 
fitting the rms of all lightcurves with low degree polynomials. 
This average rms was compared to the rms of the candidates 
with a statistical F-test, which is particularly well-suited for the 
comparison of scatter in data. All variable objects at a signifi- 
cance level of 99% are marked with a cross (and with a V in 
Table 7). 

Of the 52 CA candidates, 16 (30%) show significant vari- 
ability. In contrast, all objects in our images, i.e. mostly field 
stars, have a much smaller fraction of variable objects («10%). 
It is obvious from Fig. |5] that the fraction of variable objects 
as well as the variability amplitude does not change signifi- 
cantly over the entire magnitude range, which corresponds to 
masses down to 0.02 Mq. Figure|6]shows the same data for the 
TLS run. Unfortunately, strongly variable conditions and high 
background lead to strongly scattered rms for bright objects. 
Additionally, the high image background prevents the detec- 
tion of bright and highly variable targets, since they would be 
saturated around their lightcurve maxima. Therefore, the data 
for targets with / < 15 mag are of limited use. However, the 
plot again shows no decrease of variability amplitude with I- 
magnitude (and hence mass). Interestingly, one variable sub- 
stellar target (no. 108) was identified with an X-ray source by 
Mokler & Stelzer (2002i, confirming that this object is an ac- 
tive Brown Dwarf. 

Summarizing, these are the two results from this test: a) 
The young VLM objects in our observed field contain a large 
fraction of variable objects, b) The fraction of variable objects 
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Fig. 5. Root mean square of the time series photometry vs. I- 
band magnitude for all CA candidates. The solid line is the 
mean RMS of the relative photometry, determined by fitting the 
RMS of all objects in the field. Below 1=14 mag the long ex- 
posure photometry is influenced by saturation, brighter targets 
were analysed with the short exposure lightcurves. All objects 
marked with a cross are variable at the 99% level. 
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Fig. 6. Root mean square of the time series photometry vs. I- 
band magnitude for all TLS candidates. The solid line is the 
mean RMS of the relative photometry, determined by fitting the 
RMS of all objects in the field. All objects marked with a cross 
are variable at the 99% level. At the bright end, the variabil- 
ity test is contaminated by extinction residuals and beginning 
saturation. 



and the variability amplitude is independent of mass from 0.5 
down to 0.02 M Q . 

4.2. Period search 

Our period search is based on the widely used Scargle pe- 
riodogram (Scargle I1982> . Compared to the classical peri- 
odogram, this modified version has one advantage: The prob- 
ability that a periodogram peak at a given frequency has 
the height z should be simply P(z) — exp (—z), even for 
non-uniformly sampled datasets. This makes the Scargle pe- 



riodogram particularly well-suited for astronomical time se- 
ries with their characteristic gaps between the observing nights. 
The False Alarm Probability FAPs ca rgie for a periodic signal 
is the probability that a peak at any frequency has the height 
z: FAPscargie = 1 — [1 — exp(— z)] Ni . Ni is the number of 
independent frequencies. For evenly spaced data, Ni can be 
estimated as N 4 = -6.362 + 1.193 AT + 0.00098N 2 , where 
N is the number of datapoints (Home & Baliunas 1986 1. For 
'clumped' data, however, N (and with it the FAPs carg i e ) will 
decrease drastically. For a first estimate of the FAPscargie, we 
use Ni = N/2. Later, we will determine reliable empirical 
False Alarm Probabilities for the detected periods. 

Table 3. Candidates with significant periodic variability in 
the CA campaign. Periods (P) are determined by fitting the 
CLEANed periodogram peak with a Gaussian. Period errors 
(AP) are based on the half width at half maximum of the peri- 
odogram peak, transformed in time space. The amplitudes (A) 
correspond to the peak-to-peak-range of the binned lightcurve. 
N is the number of datapoints used for the period search. 



No. 


M(M ) 


P(h) 


AP (h) 


A (mag) 


FAP E (%) 


N 


2 


0.07 


14.7 


0.26 


0.385 


<0.01 


71 


8 


0.10 


79.1 


12.0 


0.045 


<0.01 


71 


9 


0.17 


246 


102 


0.017 


<0.01 


71 


13 


0.77 


5.63 


0.05 


0.024 


0.13 


76 


14 


0.17 


7.57 


0.09 


0.040 


<0.01 


72 


15 


0.20 


52.9 


4.24 


0.034 


<0.01 


71 


19 


0.65 


39.3 


2.34 


0.397 


<0.01 


77 


21 


0.67 


201 


74.4 


0.147 


<0.01 


76 


22 


0.06 


14.4 


0.40 


0.035 


<0.01 


69 


28 


0.03 


21.1 


0.65 


0.057 


<0.01 


71 


33 


0.17 


228 


86.4 


0.454 


<0.01 


72 


37 


0.70 


193 


66.3 


0.192 


<0.01 


78 


43 


0.06 


74.4 


8.25 


0.318 


<0.01 


72 


48 


0.11 


28.4 


1.25 


0.028 


<0.01 


71 



Table 4. Candidates with significant periodic variability in the 

TLS campaign. Columns as in Table|5] 

No. M (M e ) P (h) AP (h) A (mag) FAPe (%) N~ 

2 O07 6l6 8~14 0.547 <0.01 75" 

16 0.03 57.9 7.41 0.523 <0.01 74 

22 0.06 8.36 0.13 0.068 <0.01 76 

23 0.19 3.64 0.03 0.027 <0.01 74 

32 0.04 9.09 0.16 0.071 <0.01 76 

33 0.17 44.6 3.79 1.114 <0.01 77 
43 0.06 62.1 7.33 0.314 <0.01 77 
80 0.17 15.1 0.44 0.173 <0.01 76 
85 0.02 37.9 2.66 0.468 <0.01 74 
95 0.02 8.03 0.11 0.283 <0.01 76 

101 0.26 62.5 7.37 0.049 <0.01 75 

107 0.27 34.0 2.20 0.097 0.07 75 

128 0.02 5.78 0.07 0.066 0.39 73 



Periodograms are always contaminated with 'false' peaks, 
namely sidelobes and aliases, caused by the windowing of the 



0.07 


63.6 


0.03 


57.9 


0.06 


8.36 


0.19 


3.64 


0.04 


9.09 


0.17 


44.6 


0.06 


62.1 


0.17 


15.1 


0.02 


37.9 


0.02 


8.03 


0.26 


62.5 


0.27 


34.0 


0.02 


5.78 



8.14 


0.547 


7.41 


0.523 


0.13 


0.068 


0.03 


0.027 


0.16 


0.071 


3.79 


1.114 


7.33 


0.314 


0.44 


0.173 


2.66 


0.468 


0.11 


0.283 


7.37 


0.049 


2.20 


0.097 


0.07 


0.066 



<0.01 75 

<0.01 74 

<0.01 76 

<0.01 74 

<0.01 76 

<0.01 77 

<0.01 77 

<0.01 76 

<0.01 74 

<0.01 76 

<0.01 75 

0.07 75 

0.39 73 
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data. The recorded time series is the convolution of the sig- 
nal with this window function. Thus, any period search based 
on the raw periodogram alone will deliver spurious detections. 
To circumvent this problem, we used the CLEAN algorithm 
by Roberts et al. (19871, which deconvolves 'dirty' spectrum 
and window function and thus 'cleans' the periodogram. This 
algorithm is therefore able to distinguish between real peaks 
and spurious features in the periodogram. Since the CLEAN 
algorithm is based on the classical periodogram, it does not 
allow an a priori estimation of the FAR The combination of 
Scargle periodogram and CLEAN algorithm, however, delivers 
reliable period detections and a FAP estimate. Therefore, the 
period search procedures of numerous recent variability stud- 
ies rely exclusively on these two techniques (e.g., Bailer-Jones 
& Mundt fT999ll2U0T1 Terndrup et al. lT999l O'Dell et al. [T997l 
Patten & Simon [T9961 

As an additional test to control the significance of periods 
independently from periodograms, we propose here a method 
based on the F-test. If a detected period is not significant, then 
the variance of the time series without the period should be not 
significantly different from the variance with that period. We 
fitted each detected period with a sine wave and compared the 
scattering of the original time series a with that of the residu- 
als of the fit <7s U b. From F = cr 2 /(jg ub follows the probability 
that the two variances are equivalent. Thus, this FAPp-tost in- 
dicates, how probable it is that the found period is caused by 
chance variations in the noise of the photometry (and not in 
the periodogram). Our results (see Fig.0and|8|l show that the 
sine wave approximation is reasonable at least for low signal- 
to-noise, i.e. for objects for which reliable FAP determination 
is critical. It turns out that the FAPp-tost is in all cases larger 
than the FAP Sc argic- 

As a first step in our final period search procedure, we fil- 
tered the lightcurves to exclude obvious outliers. The following 
criteria must be fulfilled to accept a period. 

- The Scargle periodogram has a significant peak, i.e. with 
FAP Scarglc < 1%. 

- The Scargle periodograms of at least ten near neighbour 
stars show no significant peak at a similar frequency, i.e. 
no corresponding period is within ±10% of the candidate 
period. 

- Phased lightcurve and the original lightcurve itself show the 
period clearly. 

- The variance of the original lightcurve is significantly dif- 
ferent from the variance of the lightcurve without the (sine 
wave approximated) period, i.e. FAPp-test < 10%. 

- The lightcurves of nearby stars phased to the candidate pe- 
riod show no significant periodicity. 

- The CLEAN algorithm does not reject the period as spuri- 
ous feature in the periodogram. 

We established periodic variability for 14 CA and 13 TLS 
targets. Periods range from 4 to 240 hours, amplitudes from 
0.02 to 1.1 mag. Tables|3] and |4] list the relevant data for all 
objects with periodic variability. Figures[7] and [8] show the 
lightcurves phased to the measured period. 

To determine reliable FAP for the periodicities, we used 
a bootstrap approach as proposed e.g. by Kiirster et al. ([1997 1. 



The resulting empirical FAP will be called FAPe in the follow- 
ing. We obtained 10000 randomized data sets with the same 
sampling as the original lightcurve by retaining the observ- 
ing times and randomly redistributing the relative magnitudes 
amongst the observing times. We calculated the Scargle peri- 
odogram for each of these data sets and recorded the power 
of the highest peak. The FAPe is the fraction of data sets 
for which the power of the highest peak exceeds the power of 
our detected periodicity. This time-expensive simulation must 
be accomplished for every object on its own to account for 
the slightly different window function caused by our filtering. 
Therefore, we did not include these simulations in the period 
search criteria outlined above, and used instead the FAPs ca rgie 
and FAPp_tcst as preliminary significance test. It turned out, 
however, that all detected periods were indeed highly signif- 
icant with FAP E < 1%. Although the estimated FAP Scarglo 
tend to be somewhat higher than the FAPe, the two values are 
consistent within the statistical uncertainties of the simulation. 

It should be noted that the described bootstrap simulation 
possibly delivers optimistic values for the FAP. Young stars 
frequently show irregular variability on timescales of hours to 
days, and it is likely that some of our targets exhibit a similar 
behaviour. Therefore, Herbst & Wittenmyer (1996 1 discussed 
that datapoints collected within one night may not be indepen- 
dent. On the other hand, our bootstrap algorithm works on un- 
corrected data sets and makes no attempt to simulate possible 
irregular variations. Since the data sets are produced by 'shuf- 
fling' the original lightcurve, such variations are taken into ac- 
count as increased noise. 

Kearns et al. (1997 1 and Stassun et al. (1999 1 circumvent 
the problem by determining the FAP with Monte-Carlo sim- 
ulations based on synthetic lightcurves of noise with two dis- 
persions, one representing the variations during one night and 
the other representing the night-to-night variability, calculated 
as the standard deviation of nightly means. This simulation, 
however, does not adequatly reproduce intrinsic irregular vari- 
ability, because even periodic variations will cause an increased 
night-to-night scatter. Hence, it is not guaranteed that the syn- 
thetic lightcurves contain no periodic component. Therefore, 
the method will overestimate the FAP. Moreover, the method 
assumes a priori that the datapoints are normally distributed 
and that the timescale for intrinsic fluctuations is one day. A 
third drawback of the method emerges when the observing 
conditions are highly variable and there are nights with very 
few datapoints, like in our case. Under these circumstances, 
the night-to-night scatter (as defined above) is determined with 
high uncertainty. For all these reasons, we consider this ap- 
proach not applicable to our data. 

Another possibility to mimic intrinsic night-to-night vari- 
ability in simulated lightcurves is to redistribute only the inte- 
ger part of the observing dates randomly (Herbst et al. 2002 1. 
This approach, however, delivers test lightcurves with signifi- 
cantly different sampling than the original time series, a seri- 
ous drawback, because the FAP depends critically on the win- 
dow function. Moreover, this procedure only makes sense for 
time series with datapoints from many observing nights, and is 
therefore not applicable to our data. Lamm et al. (2004i com- 
pare the FAP from bootstrap simulations with and without tak- 
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Fig. 7. Phased lightcurves for the detected periodicities. No. and period from Table[3]and0]are indicated. 



ing into account the correlations between the datapoints. They 
found that the Scargle power for a FAP of 1 % decreases only 
marginally when the simulation relies on uncorrelated data- 
points. We conclude that our FAPe might be optimistic values, 
but since the FAPe for all our periods is below 0.4% (and for 
89% of our periods even below 0.01%), we are nevertheless 
confident about their significance. Finally, we note that every 
FAP can only be treated as a relative value, only comparable 
with values derived with exactly the same procedure. 



4.3. Sensitivity and completeness 

The sensitivity range of the period search is determined by our 
time sampling. For regularly spaced lightcurves, the upper fre- 
quency limit - and thus the lower period limit P m in - 1S given 
by the Nyquist frequency z^ max — j^, where A is the (con- 
stant) distance between two datapoints. The upper period limit 
Pmax, on the contrary, corresponds to the overall length of the 
time series. For irregularly spaced data, these relations are only 
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Fig. 8. Phased lightcurves for the detected periodicities (continued). 



a rough approximation. For a reliable determination of P m i n 
and Pmax, we carried out a simulation: For both campaigns, 
we selected a number of field stars which are non-variable ac- 
cording to the generic variability test described in Sect. 14.11 
Periods from 0. 1 to 300 h were added as pure sinewaves to the 
lightcurves of these objects. The amplitudes of the test periods 
were chosen so that the signal-to-noise ratio (SNR, defined as 
ratio between amplitude of periodicity and scatter in the origi- 
nal lightcurve) is larger than 2. For each period, the frequency 
of the highest peak in the Scargle periodogram was recorded. 
The difference between imposed and detected period delivers 
an estimate for the reliability of the period search for this pe- 
riod. 

Figure|9]shows the absolute difference between true and de- 
tected period for both campaigns. The dotted line corresponds 
to a period error of 10%. The simulation shows that the sen- 
sitivity varies with the period. We define P max as the longest 
period for which a detection with error below 10% is possi- 



ble and obtain 90 h for the TLS and 270 h for the CA dataset. 
Below this limit, there exist several, usually narrow windows 
where the period search with the Scargle periodogram is not 
reliable. These windows are caused by gaps in the datapoint 
distribution shown in Fig.0] Also for short periods the relia- 
bility of the period search fluctuates somewhat. This illustrates 
the commonplace that for non-continuous time series the pe- 
riod search will never be complete. In Fig.[9] we also show the 
periods found with our time series analysis (see Tables|3]andEJ 
as vertical bars. With one exception, all periods fall in regions 
where the period determination is reliable. The only exception, 
the 200 h CA period for object no. 21, will be discussed below. 

On the other hand, the simulation shows that it is easily 
possible to detect periods down to the Nyquist limit (which is 
« 0.5 h for both runs) and below. This is in agreement with 
Eyer & Bartholdi ( 1999 1 who show that for time series with 
irregular spacing the Nyquist limit is only an upper limit for 
P m i n . These results are more or less independent of the SNR, 
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Fig. 9. Sensitivity of the period search: The absolute difference 
between detected period and true period vs. true period for the 
TLS (upper panel) and the CA run (lower panel). The signal- 
to-noise ratio of the shown periodicity is 5; the dotted line cor- 
responds to a period error of 10%. The vertical bars indicate 
the periods found with our time series analysis. 



at least for values > 2, in agreement with similar simulations 
by Wolk ( 1996 1. For the simulation shown in Fig.|9] we used 
a periodicity with SNR of 5, a value typical for our periodic- 
ities. For a particular target, this plot can look slightly differ- 
ent. When the SNR of its lightcurve is very high, the width 
and heigth of the low-sensitivity windows will be slightly de- 
creased. The 200 h period for CA target no. 21 (see Fig.0 falls 
into such a window. Since the lightcurve of this object has a 
SNR of 10, much higher than the SNR of 5 for which Fig.|9] 
was computed, we consider this detection to be reliable within 
the given errors (see Table|3}. 

4.4. Pooled variance diagrams 

The period analysis described in Sect. 14. 21 detected periodic 
variability in 27 lightcurves. Some of these lightcurves, how- 
ever, clearly show variability on timescales different from the 
adopted rotation period. Hence, the corresponding phase dia- 
grams show additional variability superimposed on the adopted 
period. Moreover, some of the targets also show different pe- 



riods in both runs, as we will outline in Sect. 14. 51 To verify 
the period search results, in particular for the objects where 
the phase plots are less convincing, we therefore investigated 
the pooled-variance diagram (PVD) algorithm, which was in- 
troduced by Dobson et al. (Q990l. The basic approach of this 
method is to divide the time series into equal bins of specified 
length and to calculate the variance for each bin. The mean of 
these values (the pooled variance) measures the variability on 
the specified timescale. If the time series contains a period, the 
pooled variance should show a plateau around this period, start- 
ing roughly at P/2. Although the PVD is surely not suitable for 
precise period determination, it allows a robust verification of 
the found periods. Moreover, it will reveal if there is substantial 
evidence for variability on timescales different from the period. 
We used a PVD implementation from M. Kiirster, described in 
detail in Kiirster et al. (2000 1. The program was applied to all 
lightcurves with detected periodic variability. Especially inter- 
esting examples for the resulting PVD are shown in Fig. 1101 

The most important result of our PVD analysis is the con- 
firmation of the periods given in Tables[3] and |4] There is no 
case where the PVD shows no plateau around the adopted pe- 
riod. According to the PVD, the objects with detected periodic- 
ity clearly fall into two groups. The PVD for targets with low- 
amplitude variability reveal a broad plateau around the period, 
as shown in the upper left panel of Fig.[l0] Sometimes there is 
an additional plateau at timescales comparable with the overall 
length of our lightcurves, which can be attributed to long-term 
variations. On the other hand, objects with high photometric 
amplitude in most cases show several plateaus (see the upper 
right panel of Fig.^|for an example). Apart from a clear flat 
area around the adopted period, there are often plateaus be- 
tween 0.2 and 1 .0 days, which cannot be attributed to periodic- 
ities. Thus, for these objects the PVD indicate the existence of 
short-term variability superimposed to the period. 

The objects no. 2, 33, and 43 show high amplitude vari- 
ations in both campaigns. In the following, we discuss exem- 
plarily the PVD for object no. 2. In the lower part of Fig.^H we 
show the PVD for the CA (left panel) and the TLS (right panel) 
lightcurve. The CA lightcurve clearly shows a plateau around 
the adopted period of 14.7 h. There is also evidence for vari- 
ability on even shorter timescales, but no signs for a period of 
63.6 h, as detected in the TLS lightcurve. The PVD for the TLS 
data, however, shows several plateaus, including one around 
the 63.6 h period and a second around 16 h, in agreement with 
the period from the CA lightcurve. Hence, only the CA period 
is confirmed by both diagrams. The same result was obtained 
for objects no. 33 and 43. Since the TLS periodograms deliver 
ambiguous results, offering two possible periods, of which the 
CA periodograms confirm only one, which is also in agreement 
with the PVD analysis, we therefore consider the CA periods 
of these three targets more likely to be the rotation period. 

4.5. Long-term variability 

The variability analysis was accomplished separately for both 
campaigns. As mentioned above, 27 objects in the CA field 
were also contained in the TLS catalogue. Since both time se- 
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Fig. 10. Pooled variance diagrams for selected lightcurves. The solid lines show a running median over five datapoints. The 
periods adopted in Tables[3]and|4]are indicated as dashed line. 



ries campaigns were separated by 10 months, it is possible to 
investigate the long-term variability for those targets that were 
observed twice. The strongly different noise characteristics of 
both runs, however, complicates the comparison of the analysis 
results. 

First, we compare the generic variability test results from 
Sect. 14.11 Five objects (the high-amplitude objects no. 2, 33, 
43 and the low-level variables no. 14, 51) are variable in both 
runs. The remaining targets with variability detection in the CA 
lightcurve are either too bright or have too low amplitudes to be 
detected also in the TLS data. On the other hand, there are no 
objects which are only variable in the noisier TLS time series. 
Thus, we conclude that variability on young VLM objects can 
persist over one year. 

We examined whether a detected period is present in the 
alternative campaign as well. For nine of the doubly-detected 
targets, we find periodic variability in the CA time series. 

- The objects no. 2, 33 and 43 show high amplitude varia- 
tions in both campaigns. As discussed in Sect. 14.41 the CA 
periods can be recovered also from the TLS lightcurves, but 
not vice versa. Therefore, the CA periods are more likely to 
be the true rotation periods. 

- Targets no. 8, 9, 15, 28 show low amplitude periods in 
the CA lightcurves. None of the four objects exhibits sig- 
nificant variability in the TLS data, neither periodic nor 
non-periodic. For two of these targets, the CA amplitudes 
are too low to be detected in the more scattered TLS 



lightcurves, i.e. the average rms of the TLS lightcurves for 
the respective I-magnitude is higher than the CA amplitude. 
However, for targets 8 and 15 the CA amplitude would be 
detectable even in the TLS lightcurves. Thus, the flux mod- 
ulation must be significantly reduced compared to the CA 
data. 

- Target no. 14 has an 8 h period in the CA lightcurve and 
was significantly variable during the TLS run. However, the 
lightcurve is slightly contaminated by extinction residuals, 
possibly preventing the detection of a period. 

- Target no. 22 exhibits a significant 8.4 h period in the TLS 
lightcurve and a 14.4 h period in the CA data, i.e. the CA 
period is roughly twice the TLS period. In the next section, 
we will attribute the periodic modulations to co-rotating 
surface features. It is easily possible that the distribution of 
surface features produces a flux modulation at half of the 
rotation period, e.g. through two opposite spots. Therefore, 
the 14.4 h period is more probably the true rotation period. 



Additionally, three targets with period detection in the TLS 
run (no. 16, 23, 32) are not variable during the CA campaign. 
This is particularly surprising in the case of no. 16, since it 
shows high-amplitude variations during the TLS run. The re- 
maining two cases, no. 23 and 32, show the inverse behaviour 
of no. 8 and 15, i.e. their low-amplitude variability is only 
present in one run and at least diminished in the other. 
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5. Interpretation of the variability 

In this section, we argue that the observed variability has its ori- 
gin in two fundamentally different kinds of surface activity. We 
discriminate the two types of variability by means of the photo- 
metric amplitude, as given in Tables[3]and|4] For periodicities 
detected in the TLS and the CA lightcurve, we use the CA am- 
plitude in the following. In Fig.[^ we plot amplitude vs. angu- 
lar velocity (57 = 2ir/P) for all variable objects. Clearly, there 
is a strong cumulation of targets with low-level variability. The 
remaining objects show variability with amplitudes ranging 
from 0.28 to 0.55 mag. There is a clear gap of around 0.1 mag 
width around amplitudes of 0.25 mag, indicated by a dashed 
line in Fig.[^ Furthermore, the high-amplitude lightcurves de- 
viate visibly from a strict sine shape. To quantify these de- 
viations, we fitted all lightcurves with a sine function and 
measured the rms of the residuals. If the variability is well- 
approximated with a sine fit, the rms should not significantly 
exceed the mean rms of field stars, as given in Figs.[3]and|6] 
(solid lines). It turns out that most of the highly variable ob- 
jects with amplitudes >0.25 mag still show considerable vari- 
ability after the sine wave has been subtracted, discriminating 
them clearly from the low-amplitude objects. The classification 
in 'low-amplitude' and 'high-amplitude' variables is also con- 
firmed by the pooled variance diagrams (Sect. l4.4t , since only 
the highly variable objects show clear signs of variability on 
timescales different from the rotation period. In the following, 
we discuss possible origins for the low-level and the high-level 
variability. 
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Fig. 11. Angular velocity versus amplitude. The dashed line 
delineates the separation between low-amplitude and high- 
amplitude objects. 



5. 1. Low-level variability 

As shown above, we established periodic variability with am- 
plitudes <0.2 mag for 16 objects. Similar periodicities are well- 
known for stars in young open clusters (e.g., Rrishnamurthi et 
al. 119981 Patten & Simon [19961 , Such variability can be ex- 



plained by the existence of asymmetrically distributed surface 
features co-rotating with the objects. 

Our variable targets have effective temperatures of 
>2700K, estimated from the 3-Myr-isochrone of Baraffe et 
al. (1998 1. Hence, they are probably too hot to harbour con- 
densated dust clouds in their atmospheres (see Allard et al. 
2001 1. The origin of the low-level variability are therefore most 
probably magnetically induced spots, slightly cooler or hotter 
than their photospheric environment. The properties of these 
spots, i.e. temperature and filling factor, can only be determined 
from multi-filter/spectroscopic monitoring or future Doppler 
imaging campaigns. Nonetheless, our results give first clues 
about photospheric surface activity on very young VLM ob- 
jects. Only 17% of the targets show periodic variability - a 
significantly lower rate than in similar studies for solar mass 
stars, indicating low activity levels. This is contrary to the re- 
sults of the X-ray study by Mokler & Stelzer |2002|, who find 
no decrease of X-ray activity with mass. Their sample of X-ray 
detections could, however, be biased against objects with on- 
going accretion. Thus, at least part of the X-ray flux is possibly 
not caused by photospheric surface activity. 

There are only two objects which show low-level variabil- 
ity in both campaigns (no. 14 and 22). For these two objects 
no unambiguous period could be found. For object no. 22, this 
can be explained by evolving spot distributions (see Sect l4.5> . 
On the other hand, there are four variable objects whose flux 
modulation must be significantly reduced in one of the cam- 
paigns compared to the other, i.e. we find that the activity level 
of 50% of our objects has significantly changed in the course of 
10 months. Thus, these results give first evidence for the evo- 
lution of spot activity on young VLM objects. 

5.2. High-level variability 

Apart from the low amplitude periodic objects, there are 7 ad- 
ditional objects whose lightcurve amplitudes are >0.25 mag 
and reach up to 1.1 mag. The periodograms of their lightcurves 
show several highly significant peaks. After fitting these peri- 
ods with a sine wave and subtracting the fit, significant residu- 
als with timescales of hours remain in most cases. The ampli- 
tude of these irregular residuals is strongest for targets no. 2, 
33 and 43. We conclude that the signal of these targets consists 
of high-amplitude periodic variations, probably corresponding 
to the rotation period, and superimposed, sporadic, short-term 
features. The pooled variance diagrams for highly variable ob- 
jects imply that the timescale for these irregular variations is 
0.2... 1.0 days. (Sect.O- 

It is not possible to explain this high-level variability by 
the usual scenario of photospheric spots alone: To cause a flux 
variation of 0.75 mag, we would need a zero-emission spot with 
half the size of the visible hemisphere - this is physically im- 
possible. Eclipse events deliver no plausible interpretation as 
well, since they should cause strictly periodic light changes. 
Thus, we are confronted with a variability characteristic not 
known before for VLM objects. 

As outlined in Sect.[0 however, the observed variability 
characteristic is very similar to the high-level photometric vari- 
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ations of classical T Tauri stars (CTTS). Several groups ob- 
tained lightcurves of young solar-mass stars which show a) pe- 
riodic variability on timescales of days and b) irregular varia- 
tions on timescales of hours (e.g., Herbst et al. 2000 Fernandez 
& Eiroa 1996 Bouvier et al. 1995 1. This behaviour is usu- 
ally explained by the existence of hot spots formed by matter 
flowing from an accretion disk onto the stellar surface. Short- 
term variations can then be understood as a consequence of 
accretion rate variations and disk instabilities. Our photome- 
try shows that the photometric behaviour of CTTS extends far 
down into the substellar regime, without obvious mass depen- 
dence. This implies the existence of T Tauri-like accretion disks 
around VLM objects. 

A widely used method to check for the existence of circum- 
stellar disks is the comparison of the near-infrared colours with 
the intrinsic colours of the targets (e.g., Edwards et al. U993> . A 
disk would manifest itself by a significant near-infrared excess. 
Figure[21 shows the (H-KJ-K) colour-colour diagram for the 
cluster member candidates, constructed from 2MASS data (see 
Table 7). The solid line indicates the 3 Myr isochrone from the 
Baraffe et al. ( 1998) evolutionary tracks. Periodically variable 
objects are marked with a cross. Overplotted squares are as- 
signed to highly variable objects with amplitudes >0.25 mag, 
triangles to transition objects with amplitudes between 0.1 and 
0.2 mag. Most of the non-variable cluster members are dis- 
tributed symmetrically around the isochrone, where the scat- 
ter is determined by the photometric errors, reconfirming that 
the extinction towards a Ori is negligible and uniform over the 
cluster area (see also Bejar et al. 1999 2001 1. The periodic vari- 
able objects with low amplitude populate (with two exceptions) 
the same region in the diagram. 

The highly variable objects, in contrast, clearly tend to lie 
on the red side of the isochrone, down or even outside the red- 
dening band (dotted lines). Spectroscopy for three of these ob- 
jects assures their membership of the VLM population around 
erOri (see Sect. l7.1> . The fact, however, that the vast majority 
of the a Ori members only scatter around the isochrone shows 
that the reddened objects must suffer from intrinsic reddening. 
The fraction of objects with near-infrared excess (i.e. more than 
1 a reddening) is 7 1 % for high amplitude and 66% for transition 
objects, but 15% for low amplitude targets. Thus, we conclude 
that most of the highly variable objects exhibit a near-infrared 
excess, indicating that they are surrounded by an accretion disk. 
This may also apply to the transition objects. FigurefT2l shows 
that some accreting objects might exist among low amplitude 
and non-variable targets as well. 

The results from Sect. l4.5l show that the high-level varia- 
tions are persistent over nearly one year, although with sig- 
nificantly changed amplitude. The only exception, object no. 
16, with high-level variability in the TLS and not in the CA 
lightcurve, can be understood if it had a reduced amplitude, 
since this object has / = 18.8 mag and the mean rms at this 
brightness is 0.15 mag in the CA time series. Thus having 
shown the permanence of high-amplitude variability, these tar- 
gets are the first photometrically selected sample of strongly 
accreting VLM objects. 

Based on this conclusion, we can estimate disk frequen- 
cies and lifetimes. Out of 135 cluster member candidates, 7 
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Fig. 12. (H-KJ-K) colour-colour diagram for the cluster mem- 
ber candidates in our field. Crosses are periodically variable 
objects. Overplotted symbols flag objects with high photomet- 
ric amplitude: Squares mean that the lightcurve amplitude ex- 
ceeds 0.25 mag, and thus mark the objects defined as high- 
amplitude targets in Sect.|5] Triangles are assigned to transition 
objects with lightcurve amplitudes between 0.1 and 0.25 mag. 
The solid line indicates the 3 Myr isochrone from the Baraffe 
et al. evolutionary tracks. Dotted lines show the interstellar ex- 
tinction vector calculated from Mathis (1990). A reddening 
vector for Ay = 1 mag is given. For some reddened objects, 
we overplot the errorbars for the photometry. 

show high-amplitude variations and thus probably possess an 
accretion disk. This corresponds to a fraction of 5 to 7% (tak- 
ing into account a field star contamination of up to 30% of our 
VLM sample). Since we might have missed accreting objects 
without strong photometric variability, this should be consid- 
ered as a lower limit. For comparison, Barrado y Navascues et 
al. ( I2003i . Oliveira et al. d200lt . and Barrado y Navascues & 
Martin (2003 1 estimate disk frequencies of 5-14% for VLM 
objects in the a Ori cluster. Based on a very small sample of 
six Brown Dwarfs, Jayawardhana et al. (2003 1 derived a disk 
frequency of 33 ± 24%. Given the statistical uncertainties, all 
these results are in good agreement with our estimate. Thus, it 
seems that disks are rare in the a Ori cluster, but they exist for 
VLM stars as well as for Brown Dwarfs. The fact that we do 
not observe totally irregular variability with high amplitudes is 
a further indication for the scarcity of strong accretors among 
VLM objects in the a Ori cluster, since such variations are be- 
lieved to be typical for the most actively accreting T Tauri stars 
in younger clusters (e.g., Herbst et al. 2000). 

The low disk frequency for VLM objects in the a Ori clus- 
ter implies that most of these objects dissipate their disks be- 
fore they reach the age of this cluster. Since the most probable 
age for the cluster is 3 Myr (Zapatero Osorio et al. 2002 1, the 
disk lifetime for VLM objects, defined as the timescale for es- 
sentially all the stars to lose their disk (Haisch et al. 2001}, is 
roughly 3-4 Myr. In comparison, solar-mass stars retain their 
disks significantly longer (Haisch et al. 120011 . their disk life- 
time is 6 Myr. Thus, if the age estimate for the a Ori cluster 
is valid, disk dissipation timescales are shorter in the VLM 
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regime. This is confirmed by a recent study of Oliveira et al. 
( 2003 1, who investigated near-infrared colours of a sample of 
crOri cluster members with J < 13.3 mag (corresponding to 
M > 0.2 Mq). They derived an overall disk frequency of 46- 
54%, but the near-infrared excess shows a strong decline for 
K > 11 mag, corresponding to M < 0.5 Mq (see their Fig. 
5). Thus, most of the stars with disks have solar-like masses, 
and the disk frequency is decreasing with mass. 

6. Rotation Periods 

As discussed above, we attribute the detected periodicities to 
a rotational modulation of the emitted flux. In this section, we 
will analyse the implied rotation periods. We will compare the 
periods from high-level variable with those for low-level vari- 
able objects. We will then discuss the mass dependence of rota- 
tion. The period sample for all analyses comprises 14 CA and 
9 TLS periods, excluding those TLS periodicities also detected 
in the CA lightcurves (see also the discussion of these objects 
in Sect. l4.5> . The two rotation periods for a Ori objects (SOri 
31 and SOri 33) measured by Bailer-Jones & Mundt J200H 
will be used to underline the results. 

6. 1. High amplitude vs. low amplitude 

In Sect.|5]we argue that we observe two kinds of variability, 
high-level and low-level variations, where the high-level vari- 
ability is interpreted as a consequence of ongoing accretion 
processes, as confirmed by near-infrared photometry (Sect. l5.2l 
and spectroscopy (see Sect.0. Thus, highly variable objects 
are very likely surrounded by circumstellar disks. This could 
have consequences for their rotational behaviour, since mag- 
netic coupling between star and circumstellar disk is believed 
to be the mechanism responsible to brake the rotation of T Tauri 
stars. For solar mass stars, there is indeed observational evi- 
dence for a correlation between rotation and disk indicators, 
in the sense that stars with disks tend to be slow rotators (e.g. 
Edwards et al. fl993l Herbst et al. I20021 . However, there ex- 
ist controversial investigations which find no correlation be- 
tween rotation and disk presence (e.g., Stassun et al. 20THV in- 
dicating that disk-locking may not be the full solution to angu- 
lar momentum distribution in young stars (Rebull et al. 2002 1. 
Nevertheless, assuming that disk-locking plays a major role as 
an angular momentum regulation mechanism, we should ex- 
pect that accreting VLM objects show a rotation period dis- 
tribution different from the non-accreting ones. In particular, 
objects with disks should on average rotate more slowly. 

In Fig. ^3 we plot photometric amplitude vs. angular ve- 
locity for all our targets. As argued in Sect. l5.2l the plot shows 
two populations of datapoints, separated by a gap around am- 
plitudes of 0.2 . . . 0.25 mag, delineated by a dashed line. The 
objects above this line are probably active accretors. For these 
objects, the average angular velocity is ft = 6.0 dr 1 . For am- 
plitudes < 0.25 mag, however, ft is on average 11.1 dr 1 . Thus, 
low- amplitude objects tend to rotate faster on average, as we 
would expect for a disk-locking scenario. The periods from 
Bailer- Jones & Mundt (2001) are also consistent with this pic- 
ture, since they are shorter than 10 h (corresponding to O > 15) 



and have amplitudes < 0.05 mag. Hence, the available data 
may indicate disk-locking processes on VLM objects. Since 
this result could be affected by low number statistics, we post- 
pone its further discussion until more rotational data are on 
hand. 

It is instructive to compare our rotation periods with those 
recently published by Joergens et al. (2003 1. They measured 4 
periods for Brown Dwarfs and VLM stars in the Chal star form- 
ing region, which is 1 Myr old and thus considerable younger 
than the crOri cluster. All their periods are longer than 2 days. 
In the context of the disk-locking scenario, these objects could 
be still disk-locked. This would be in agreement with the large 
disk frequencies (65%) among Brown Dwarfs in the similarly- 
aged Trapezium cluster (Muench et al. l2001i . 

6.2. Mass dependence 

Several groups have reported a bimodal distribution of the ro- 
tation periods for very young solar mass stars (e.g. Herbst et 
al. 2002 Edwards et al. 1993 1. Since longer periods are signifi- 
cantly correlated with infrared excess attributed to the presence 
of a circumstellar disk, this is a direct sign for a disk-locking 
mechanism braking the rotation in the T Tauri phase. There is 
observational evidence, however, that the bimodality becomes 
less obvious for low-mass objects and finally vanishes for ob- 
jects with masses below 0.2 Mq (Herbst et al. 2001 1. Our pe- 
riods show no bimodal distribution throughout the whole mass 
regime. Since we observed only few targets with masses com- 
parable to the solar mass, the detection of the slow rotator peak 
in the histogram could be prevented by our low-number statis- 
tics. Hence, our data allow no verification of the evolution of 
bimodality with decreasing mass. Nevertheless, we note that 
there is definitely no bimodality within our statistical uncertain- 
ties for objects with masses in the range from 0.03 to 0.4 Mq. 

To quantify the mass dependence of the periods, Herbst et 
al. (|2001 1 calculate a moving median for their periods in the 
mass regime 0.1 .. . 0.4 Mq, measured for objects in the ONC. 
They find that lower-mass objects tend to rotate gradually 
faster. Our results confirm this result. FigurelT3lshows masses 
vs. rotation periods for our a Ori targets. The median period is 
46.1 h for VLM stars, but only 14.7 h for Brown Dwarfs. Our 
periods make it possible for the first time to extend the analysis 
of Herbst et al. (2001 1 to masses below their detection limit. 
Figure[T3l shows the mass-period relation of Herbst et al. com- 
bined with our datapoint for masses < 0.1 Mq. Obviously, our 
periods nicely extend the mass-rotation relation down into the 
substellar regime. Because of the statistical uncertainties, how- 
ever, this needs further confirmation. 

Recapitulating, we observe two effects: no bimodality and 
the tendency of faster rotation towards lower masses. Both 
could be explained by a decrease of the effectiveness of ro- 
tational braking through magnetic coupling to a circumstellar 
disk. Indeed, there is evidence for decreasing disk lifetimes 
with decreasing mass, as pointed out in Sect. 16.11 Moreover, 
efficient angular momentum removal, either through magnetic 
coupling or stellar winds, could also be prevented by low mag- 
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Fig. 13. Mass dependence of the rotation period: Plotted are 
all measured periods from TLS and CA vs. the object mass. 
The solid line shows the median period published by Herbst 
et al. 



J2001 1. The median for our objects with M < 0.1 M© 
is shown as filled square, vertical lines indicate the quartiles 
(giving an estimate of the statistical uncertainties). This median 
found for our objects extends the line of Herbst et al. (|2001 ) 
into the Brown Dwarf regime (dashed line). 

netic field strengths, as implied by the the low activity level 
inferred from the small rate of variable objects (see Sect. 15. II . 

7. Accreting VLM objects 

In this section, we concentrate again on the discussion of the 
highly variable targets, which were identified as actively ac- 
creting VLM objects in Sect.[3] To substantiate this interpre- 
tation, we obtained low-resolution spectroscopy with MOSCA 
at the 3.5-m telescope on Calar Alto and grism red500, cover- 
ing a wavelength range of 5400 . . . 10000 A with a dispersion 
of 2.9 A/pixel. We observed two Brown Dwarfs and one VLM 
star with high-level variability in both time series (no. 2, 33, 43) 
and, for comparison, one non-variable Brown Dwarf (no. 90). 
For relative flux calibration, we additionally obtained a spec- 
trum of the spectrophotometric standard star GD71 (Bohlin 
et al. I1995i . Spectra of Ar/ArHg/Ne arc lamps were used for 
wavelength calibration. Reduction and calibration was done 
with the IRAF standard routines within the package onedspec. 
We note that the flux is considerably damped for A > 8000 A, 
because of second order overlap in the spectra. All four spectra 
are shown in Fig.1141 

7.1. Spectral classification 

The spectra show all the characteristics expected from late-type 
objects, in particular deep TiO and VO bands (Kirkpatrick et 
al. Q991). We determined T c g and \ogg from the spectra by 
comparing them with the STARdusty2000 models of Allard 
et al. J2000> . For this comparison, we used only the spectral 
range from 6700 to 7700 A, to exclude regions with strong 
emission lines and avoid the long wavelength end, where our 
spectra are not properly flux calibrated. The model grid spans 
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Fig. 14. Optical spectra for three highly variable targets (from 
top to bottom: no. 2, 33, 43) and one non-variable target (no. 
90, beneath). The spectra are shifted by 12, 6, 4, (from top to 
bottom) units for clarity. 

T cff = 2000... 4000 Kin 100 K steps and logg = 3.5... 5. 5 
in steps of 0.5. To account for varying atmospheric conditions 
during the observations and for possible continuum excesses as 
expected for accreting objects, we allowed for a wavelength- 
independent scaling and shifting of the models. 

We calculated the root mean square between each model 
and the data. It turns out that the rms values strongly depend on 
T e g, with a clear minimum. On the other hand, the dependency 
on log g is much weaker. Therefore, we first determined the 
effective temperature at which the rms is smallest. From the 
set of models with this T e g, we then selected the model with 
the best fitting log g. Visual comparison of the selected model 
with the data shows good agreement. The results, i.e. T e g and 
logg, are listed in Table|5] Spectral types are estimated with 
the effective temperature scale of Leggett et al. (2000 1, which 
is based on essentially the same models as our analysis. 

Table 5. Comparison of the spectra with models of Allard et al. 
( 2000). The first two columns are the parameters of the best fit- 
ting model spectrum. The spectral types follow with the effec- 
tive temperature scale of Leggett et al. ( 2000!. F° r comparison, 
we additionally give the T c g estimated from photometry alone, 
presuming cluster membership. * Non-variable target observed 
for comparison 



no. 


2 


33 


43 


90* 


T eff (K) 


2800 


3100 


2900 


3000 




4.0 


3.5 


4.0 


3.5 


Spectral type 


M5.5 


M4.0 


M5.0 


M4.5 


T cff (K)-phot 


2930 


3030 


2890 


2850 



Assuming that the object is indeed a cluster member, a sec- 
ond estimate for T c g can be obtained by comparing the pho- 
tometry with the cluster isochrone of Baraffe et al. ( 1998 1. This 
value of X" ff will be either too high or too low, if the target 
resides in the fore- or background of the cluster. An estimate 
in agreement with the true T c g, as measured from the spec- 
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tra, confirms the cluster membership. We cannot expect perfect 
agreement, however, since this comparison is affected by pho- 
tometry errors and the grid size of the models. A further com- 
plication is of course the large variability of targets no. 2, 33 
and 43. It might be that photometry and spectroscopy caught 
the object in a totally different accretion state. With these con- 
siderations in mind, all objects show good agreement between 
spectroscopic and photometric effective temperatures, confirm- 
ing them als VLM cluster members of crOri. The youth of all 
four targets its also confirmed by the low surface gravity values. 
Evolved objects with comparable masses exhibit \ogg > 5.3 
(Baraffe et al. 1998 1, whereas we derived \ogg = 3.5 . . . 4.0, 
in agreement with the values expected for a 3 Myr object with 
T eS around 3000 K. 

7.2. Emission-line spectrum 

The dominant spectral feature in all three highly variable ob- 
jects is the large Ha emission line. The non-variable object has 
Ha in emission as well, but with considerably lower intensity. 
A variety of other emission features is present in the spectra, 
similar to emission spectra of CTTS. We measured equivalent 
widths for all prominent emission lines by subtracting a lin- 
early fitted continuum from the line flux (Table|6j. Obviously, 
the highly variable targets, in particular targets no. 2 and 43, 
show much more intense emission features than the compari- 
son object. Target no. 33 seems to be more active than no. 90, 
judged from the Ha equivalent widths. However, compared to 
no. 2 and 43, Ha and Ca emission lines are relatively faint. It 
might be that we caught target no. 33 at the minimum of its 
lightcurve and thus in a relatively quiet state. 

Forbidden lines such as [OI] and [SII] can be attributed to 
low-density regions, e.g. in stellar winds, and are thus tracers 
of a mass loss process. On the other hand, features like the Ha 
emission line and the Ca triplet are clear signs of the accretion 
process. These lines are also common in outflows. For both fea- 
tures, however, we attribute most of the emission to accretion 
processes: On the one hand, the ratio between Ha and [SII] 
in our spectra is much larger than 4, which is the upper limit 
given for (outflow-dominated) Herbig Haro objects (e.g. Bohm 
& Goodson 1997 1. On the other hand, the ratio of the Ca triplet 
equivalent widths is close to 1:1:1, whereas the predicted ra- 
tio for optically thin emission would be 1:9:5. Similar values 
are common among strongly accreting CTTS (e.g., Reipurth et 
al.[H)86l and clearly differ from those of Herbig Haro flows. 
The presence of the HeA6678 line is a third indicator of ongo- 
ing accretion. We conclude that the emission spectrum of those 
highly variable objects is dominated by the signature of strong 
accretion processes. Comparing Ha linewidths and spectral 
types with the classification scheme of White & Basri (2003 1, 
the highly variable objects belong to the VLM analogues of 
classical T Tauri stars, while the non-variable object is similar 
to a weak-line T Tauri star. 

Recapitulating, our low-resolution spectroscopy yields a 
clear result: Our highly variable VLM objects are probably ac- 
creting material from a disk and are thus the VLM analogues 
of CTTS. The T Tauri nature of young VLM objects is sup- 



Table 6. Equivalent widths for prominent emission lines. The 
first column contains name and laboratory wavelength (in A) 
of the emission line, then we list the equivalent widths for the 
four targets. Non-detection of the feature is indicated with a 
dash, 'b' means that the line is present, but strongly blended by 
nearby features. 

* Non-variable target observed for comparison 



no. 


2 


33 


43 


90* 


[OI] 6300 


b 


1.1 


b 




[OI] 6364 


0.9 








Ha 6563 


104.5 


21.3 


103.7 


10.7 


Hel 6678 


1.8 


2.4 






[SII] 6716 


0.23 


0.39 






[SII] 6731 


0.38 


0.32 






OI 8446 


b 








Call 8498 


4.8 


0.6 


6.9 


3.8 


Call 8542 


3.8 




5.1 




Call 8662 


2.5 




3.1 




HI (P9) 9229 


0.8 


0.4 


0.2 


0.4 



ported by recent publications, since spectroscopic signatures of 
typical T Tauri behaviour have been detected by other groups 
(e.g., Jayawardha na et a l. 120021 White & Basri 2003 Barrado 
y Navascues et al. 2003 ). The existence of disks around Brown 
Dwarfs is confirmed by independent observations of Natta & 
Testi J200TI . Nat ta et a l. d20U2l . Testi et al. (120021 , Apai et al. 
( 2002|i, Liu et al. J2003t . Jayawardhana et al. J2003> . and Lopez 
Marti et al. (2004i, mostly via mid-infrared excesses. All these 
findings set strong constraints for possible formation scenar- 
ios for VLM objects. In competition with star-like formation 
mechanisms, it was recently proposed that Brown Dwarfs form 
preferably in orbits around more massive stars, either in pro- 
toplanetary disks (Pickett et al. 2000|l or as stellar embryos in 
multiple stellar systems (Reipurth & Clarke 2001). Both sce- 
narios need an ejection mechanism to explain the large num- 
ber of isolated Brown Dwarfs found in stellar clusters and star 
forming regions. It seems unlikely that the objects maintain a 
significant fraction of their accretion reservoir during this ejec- 
tion process. Thus, the current observational picture favours 
star-like formation for VLM stars and substellar objects. 

8. Conclusions 

We present the first photometric variability study for VLM 
members of the young a Ori cluster. With multi-filter photom- 
etry, 135 cluster member candidates were identified, including 
90% VLM objects (M < 0.4 M©). This preliminary member 
list still contains about 30% contaminating field stars and thus 
needs spectroscopic verification. 

We monitored these targets in two I-band time series in 
January and December 2001 . The first campaign, using the 2-m 
Tautenburg Schmidt telescope, covered 110 candidates, span- 
ning masses from 0.02 to 0.5 M@. The second time series, ob- 
tained with the 1.23-m telescope on Calar Alto, registered 52 
targets spanning a mass range from 0.02 to 1 .4 Mq . 27 Objects 
were observed in both runs. After reduction and relative cali- 
bration, we obtained differential lightcurves for all targets. We 
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found that the difference image analysis package from Gossl et 
al. (|2002 1 improves the precision of the Calar Alto lightcurves 
by several mmag. 

Time series analysis was focused on the period search, but 
we included a test to detect irregular variations as well. We 
examined the variability of the targets based on the scattering 
of the lightcurves and found that the very young targets contain 
significantly more variable objects than field stars. The fraction 
of variable objects as well as the variability amplitude does not 
significantly change over the entire mass range. 

We detected periodic variability for 14 CA and 13 TLS tar- 
gets, with periods spanning from 4 to 240 hours. Four targets 
show a periodicity in both campaigns. The periodically vari- 
able candidates clearly fall into two groups: 16 of them show 
low-level variability which is in most cases well-approximated 
by a sine wave. We argue that this type of variability is most 
probably caused by magnetically induced photospheric spots 
co-rotating with the objects. We see evidence for the spot evo- 
lution on VLM objects, since most of these targets show the 
periodicity in only one campaign. 

The lightcurves of the remaining objects have amplitudes 
>0.25 mag and show obvious deviations from the sine shape. 
Their photometric behaviour is very similar to those of classical 
T Tauri stars. The quasi-periodic oscillation is therefore most 
probably caused by hot spots formed by matter flow from an 
accretion disk. Thus, our results imply the existence of accre- 
tion disks around VLM stars and Brown Dwarfs. The fact that 
most of the highly variable objects show a near-infrared colour 
excess confirms this finding. Low-resolution spectroscopy of 
a subsample of these targets assures their cluster membership 
and their youth. We find strong emission features, in particu- 
lar a dominating Ha emission line, characteristic for accreting 
objects. The extension of T Tauri behaviour to VLM stars and 
Brown Dwarfs strongly supports star-like formation scenarios 
for these objects. 

Highly variable and thus accreting objects rotate on aver- 
age more slowly than low-level variable targets. This is ex- 
pected in terms of a disk-locking scenario, where magnetic 
coupling between star and disk removes angular momentum 
during the T Tauri phase. In agreement with previous publica- 
tions, we find that strong accretors are rare in the a Ori clus- 
ter, but they certainly exist among VLM stars as well as Brown 
Dwarfs. Comparison with studies of younger clusters and more 
massive stars give evidence that disk dissipation timescales are 
significantly decreased with decreasing object mass. 

In agreement with previous publications, we find a ten- 
dency of faster rotation towards decreasing object masses, pos- 
sibly caused by decreasing effectiveness of angular momentum 
removal. This could either be explained with shorter disk life- 
times or with weaker magnetic fields, as inferred from the low 
activity level of our targets. 
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